Mirrorless Negative-index Parametric Micro-oscillator 
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The feasibility and extraordinary properties of mirrorless parametric oscillations in strongly ab- 
sorbing negative-index metamaterials are shown. They stem from the backwardness of electromag- 
netic waves inherent to this type of metamaterials. 
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Optical negative-index (NI) metamaterials (NIMs) 
form a novel class of artificial electromagnetic media 
that promises revolutionary breakthroughs in photon- 
ics (see for example [l|). Unlike ordinary positive- index 
(PI) materials (PIMs), the energy flow and wave vec- 
tor are counter-directed in NIMs, which determines their 
extraordinary linear and nonlinear optical (NLO) prop- 
erties. Strong absorption is fundamentally inherent to 
NIMs, which presents a severe detrimental factor toward 
their applications. Herein, we show the feasibility and ex- 
traordinary properties of mirrorless generation of entan- 
gled counter-propagating left- and right-handed photons 
in a strongly absorbing micro sample of a NIM. 

Frequency-degenerate multi-wave mixing and self- 
oscillations of counter- propagating waves in ordinary ma- 
terials have been extensively studied because of their 
easily achievable phase matching. Phase matching for 
three- wave mixing (TWM) and four- wave mixing (FWM) 
of contra-propagating waves that are far from degener- 
acy seem impossible in ordinary materials. The possi- 
bility of TWM with mirrorless self-oscillations from two 
co-propagating waves with nearly degenerate frequen- 
cies that fall within an anomalous dispersion frequency 
domain and a far-infrared difference-frequency counter- 
propagating wave in an anisotropic crystal was :>posed 
in [2] (and references therein) more than 40 ,rs ago 
(see also |3|). However, far- infrared radiation is pically 
strongly absorbed in crystals. For the first tin TWM 
backward-wave (BW) mirrorless optical param cal os- 
cillator (BWMOPO) with all three different opt [ wave- 
lengths was recently realized [4|. Phase-ma ing of 
counter-propagating waves has been achieved in a sub- 
micrometer periodically poled NLO crystal, which has 
become possible owing to recent advances in nanotech- 
nology. Both in the proposal [2j and in the experiment 
[4| , the opposite orientation of wave vectors was required 
for mirrorless oscillations due to the fact that PI crys- 
tals were implemented. As outlined, a major technical 
problem in creating BWMOPO stems from the require- 
ment of phase matching with the opposite orientation of 
wave vectors in PIMs. Herein, we show the feasibility 
of creating distributed feedback and BWMOPO where 



an antiparallel orientation of wave vectors of the coupled 
waves is not required anymore. Thus BWMOPO at ap- 
preciably different frequencies becomes possible while all 
the wave vectors of the coupled waves remain co-directed. 
Such an opportunity makes phase matching much easier 
and is offered by the backwardness of electromagnetic 
waves, which is natural to NIMs. We note that NLO in 
NIMs still remains a less developed area of electromag- 
netism. Two options are proposed in this work. One is 
TWM BWMOPO, which implements intrinsic x^ non- 
linearities associated, for example, with the asymmetry of 
current-voltage characteristics from the building blocks 
of NIMs [5j • A recent demonstration of the exciting pos- 
sibilities to craft NIMs with NLO responses that exceed 
those from natural crystals is reported in [6| . This option 
allows broad frequency tunability; however, there might 
be technical problems to overlap the NI, phase match- 
ing and strong NLO response frequency domains. The 
second option, FWM BWMOPO, offers independent en- 
gineering of a strong nonlinearity x through embed- 
ded resonant NLO centers. In the vicinity of resonances, 
it becomes exceptionally strong. In addition, the opti- 
cal properties of the composite become tailored by the 
means of quantum control. 
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FIG. 1: Coupling geometry for three-wave BWMOPO, (a); 
the tailored distribution of the signal, T(z) = \hi(z)/hii,\ , 
and of the idler, r/2 — \E2 / Ei{L)\ 2 , along the slab at Afc = 0, 
(b) and (c). (b) and main plot (c): a±L — 2.5, aiL = —2. (b): 
main plot: gL = 7.865; upper right inset: four-order decrease 
in Ti(z) outside the geometrical resonance at gL = 10; upper 
left inset: geometrical resonances in output signal T\(z — 0). 
(c): gL — 1.554. Upper left insert: 01.2L = — aiL = —2.5. 
Upper right insert: ct^L = a±L = 2.5. 



Figure HJa) depicts the coupling geometry for the pro- 
posed TWM BWMOPO. We assume the wave at u x with 
the wavevector ki directed along the z-axis to be a NI 
(rii < 0) signal. Therefore, it is a backward wave be- 
cause its energy-flow Si = (c/47r)[Ei x Hi] appears di- 
rected against the z-axis. The medium is illuminated by 
a higher-frequency PI wave at U3 traveling along the pr- 
axis (n.3 > 0). The two coupled waves with co-directed 
wave vectors k3 and ki generate a difference-frequency 
idler at 0J2 = W3 — Wi j which is also assumed to be a PI 
wave (ri2 > 0). The idler contributes back into the wave 
at <ji\ through TWM and thus enables optical paramet- 
ric amplification (OPA) at uj\ by converting the energy 
of the control field at W3 into the signal. Thus, all of the 
three coupled waves have their wave vectors co-directed 
along z, whereas the energy flow of the signal wave, Si, 
is counter- directed to the energy flows of the two other 
waves, S2 and S3. Such a coupling scheme is in strict con- 
trast both with the conventional phase-matching scheme 
for OPA in ordinary crystals, where all energy-flows and 
phase velocities are co-directed, and with the TWM BW- 
MOPO outlined above, where the energy flow and wave 
vector of one of the waves are opposite to all others. The 
equations for slowly- varying amplitudes of the signal and 
the idler take the form 

dhi/dz = iaih^h^ exp[zAfcz] + (a.\/2)h\, (1) 

dhijdz = io~2hsh\ exp[iAfcz] — (a2/2)/i2, (2) 

Here, hj are the magnetic components of the fields, kj — 



\rij\u)j/c > 0; <7j = •iTTx^Sj^'j/kjC 2 ; Ak = fc 3 — fc 2 



fci; 



X~ ' is the magnetic nonlinear susceptibility; Oj are ab- 
sorption indices. The magnitude /13 is assumed con- 
stant along the medium. Equation (JlJ exhibits three 
fundamental differences as compared with TWM of co- 
propagating waves in ordinary materials. First, there is 
an opposite sign of <j\ because of t\ < 0. Second, there is 
an opposite sign before ct\ because the energy-flow Si is 
directed against the axis z. Third, the boundary condi- 
tions for hi are defined at the opposite edge of the sample 
as compared to hi and /13 because the energy-flows Si 
and S2 are counter-directed. The transmission factor for 
the negative-index signal, T(z) = \hi(z)/hn,\ , is derived 
from the solution to the equations (JTJ) and |(2J| as [7( 



Ti(z = 0) = T w = 



exp{-[(ai/2)- S ]L} 



cos RL + (s/R) sin RL 



(3) 



Here, s = (en + a 2 )/(4) - i(Afc/2), R = y/g 2 - s 2 and 

g = ( % /wiw 2 /\/eie2/MiA f 2)(87r/c)x' 2 ^3- The fundamen- 
tal difference between the spatial distribution of the sig- 
nal in ordinary and NI materials is explicitly seen at 
a. 



x 3 — Ak = 0. Then equation ([3|) reduces to 
T w = l/[cos(gL)} 2 , 



(4) 



whereas, in ordinary media, the signal would exponen- 
tially grow as Ti ex exp(2gL). Equations (O and Q 



present a sequence of geometrical resonances as func- 
tions of the slab thickness and of the intensity of the 
control field /13. Such extraordinary resonances provide 
for the feasibility of attaining the oscillation threshold 
for the generation of the entangled counter-propagating 
left-handed, Hu>i, and right-handed, h\oi, photons with- 
out a cavity at gL — > (2j + l)ir/2. A similar behavior is 
characteristic for distributed-feedback lasers and is equiv- 
alent to a great extension of the NLO coupling length. 
It is known that even weak amplification per unit length 
may lead to lasing provided that the corresponding fre- 
quency coincides with a high-quality cavity or feedback 
resonances. The crucial role of the outlined geometrical 
resonances, as well as the dramatic change caused by am- 
plification of the idler, is illustrated in Fig. Up)), (c). The 
main plot in Fig. [IJb) shows the spatial behavior of the 
signal and the idler at a slightly off-resonant value of gL, 
which results in the overall decrease of the transmitted 
signal. Since the idler grows toward the back facet of the 
slab at z = L and the signal experiences absorption in 
the opposite direction, the maximum of the signal for the 
given parameters appears closer to the back facet of the 
slab. In general, the distribution of the signal along the 
slab and, consequently, its output value depends strongly 
on the difference of the absorption indices for the signal 
and the idler. As seen from the comparison with the up- 
per right inset, the signal maximum and the slab trans- 
parency drop sharply with offset (here, an increase) in 
the magnitude of gL from its resonant value. The upper 
left inset in Fig. [ljb) depicts the resonances in the out- 
put signal at z = 0. A change in the slab thickness or 
in the intensity of the control fields leads to significant 
changes in the distributions of the signal and idler along 
the slab, which is seen from a comparison of the main 
plots in Fig. QJb) and (c). Such a dependence is in stark 
contrast with its counterpart in PI materials. The role of 
the additional amplification of the idler provided by the 
control fields is readily seen from the comparison of the 
upper left and right insets in Fig. QJc) . 

The basic idea of resonant FWM BWMOPO is de- 
picted in Fig. [2fa),(b). A slab of NIM is doped by four- 
level nonlinear centers [Fig. Ufa)] so that the signal fre- 
quency, ll>i, falls in the NI domain (n(u>i) < 0), whereas 
all the other frequencies, 0J3, LU4 and u)%, are in the PI 
domain. Here, we look for the possibility to achieve the 
oscillation threshold for the BW at wi controlled by two 
lasers at lo^ and W4. Due to FWM, these three fields 
generate an idler at 0J2 = ^3 + <^4 — wi , which may expe- 
rience significant amplification, either Raman or due to 
the population inversion at the corresponding resonant 
transition, driven by the control fields. This opens an 
additional channel that may greatly enhance the coher- 
ent energy transfer from the control fields to the signal. 
The amplified idler contributes back to Wi = lo^, +LU4 — u>2 
through FWM and thus causes strongly enhanced para- 
metric amplification of the signal. Unlike ordinary off- 




FIG. 2: Scheme of quantum-controlled resonant FWM BW- 
MOPO, (a); corresponding coupling geometry (m < 0), (b); 
laser-induced nonlinear interference structures and dispersion 
of local parameters, (c)-(e); transmittance T\(z = 0) of the 
NIM slab vs Ke(gL) and lm(gL) at different phase mismatch 
8kL, (f) and (g); BWMOPO threshold vs frequency resonance 
offset and the NIM slab thickness, (h). (c)-(e): Gz = Ga = 
20 GHz, ?L 3 = -fi 4 = T gt . y = Q.i/Ti n . Solid: Im( 5 /ai ) 2 ; 
dash: Re(g/aio) 2 ; dash-dot: [(02 + ai + ajv/A/)/4Qio] 2 . (c): 
overview, (d): left resonance, (e): right resonance. a^L = —2, 
a x L = 2.5. (f): Afc = 0. (g): AkL = -7tt. (h): 
Sk/a 10 = 0.248, L r = 32.6426L™, y =-53.13. 



resonant NLO, in this scenario many-orders of resonance 
enhancement in the NLO coupling is accompanied by a 
strong change of the local optical parameters by the con- 
trol fields. This is because the control fields may pro- 
duce a significant population transfer and even inversion, 
as well as the modulation of the probability amplitudes. 
Consequently, a split and other modifications of the reso- 
nance shapes occur that stem from the quantum interfer- 
ence, which is shown in Fig. [DJc)-(e). Switching between 
constructive and destructive interference may serve as 
a tool for harnessing local optical coefficients. Alterna- 
tively, such changes can be minimized, so that the major 
amplification would come directly from the coherent en- 
ergy transfer from the control fields to the signal and the 
idler through the FWM coupling. The interference of 
quantum pathways in the vicinity of the resonances may 
even lead to the fact that the overall process ceases to be 
seen as a set of successive one- and multi-photon elemen- 
tary processes [9j (and references therein). Similar to the 
case for TWM, the features of FWM in NIMs employed 
here appear in stark contrast with those known for the 
ordinary PI materials due to the fact that all wave vectors 
are co-directed in order to ensure maximum phase match- 



ing but the energy flow Si appears contra-directed to I44 
and, therefore to all other wave vectors [Fig. |5fb)]. Like 
in the preceding case, this imposes extraordinary features 
on the nonlinear propagation of the coupled waves. The 
slowly-varying electric amplitudes of the coupled waves 
at u>i and u>2 are given by the equations 

dE 1 /dz ~ -i 1 [ 3) E'*exp[iAkz] + (a 1 /2)E 4 , (5) 

(6) 



dE 2 /dz = i^ 3) E$ exp[iAkz] - (a 2 /2)E 2 . 



Here, 7JJ = (47r|/ii,2|wi 2 /ki,2C 2 )x[ 2E3E4 are NLO cou- 
pling coefficients; the dielectric permittivity and mag- 
netic permeability, ej and jij, are negative at w\; Afc = 
&3 + ki — k\ — ^2', cij are the absorption or amplifica- 
tion coefficients; lj\ + u>2 = ^3 + ^4; kj = \n,j\ujj/c > 0; 

(3) 
and x\ 2 are effective electric nonlinear susceptibilities. 

The amplitudes of the fundamental (control) waves £3 
and E4 are assumed constant along the slab. Equations 
§5$l and ^ are similar to equations ([T]) and ^ and ex- 
hibit the same fundamental differences as compared with 

their counterparts in ordinary PI materials: the sign of 

(3) 
the nonlinear polarization term 7J is opposite to that 

(3) 

of 72 , which occurs because fii < 0; the opposite sign 
appears for ai; and the boundary conditions for E\ must 
be defined at the opposite side of the sample as com- 
pared to those for all other waves. These differences 
lead to a counterintuitive evolution of the signal and idler 
along the medium, similar to those discussed for the case 
of TWM BWMOPO. With the aid of the solutions to 
equations ^ and © [8], the transmission (amplifica- 
tion) factor for the NI signal, T x {z) = \E 1 (z)/E 1L \ 2 , is 
given by equation l[5]). where g 2 = 72 7i- A significant 
difference between the resonant and off-resonant NLO 
processes investigated here is that the NLO susceptibil- 
ities and, therefore, the parameters 71 and 72 become 
complex and different from each other in the vicinity of 
the resonances. Hence, the factor g 2 may become nega- 
tive or complex. This indicates an additional phase shift 
that causes further radical changes in the nonlinear prop- 
agation features which can be tailored. Figures H{c)- 
(e) show such laser-induced nonlinear interference res- 
onances in the coupling parameter g 2 with commensu- 
rate real and imaginary parts. Figures EJf),(g) prove 
that resonance coupling and, hence, complex or nega- 
tive magnitudes of g 2 may result in the fact that phase 
matching Afc = ceases to be required, and the large 
phase mismatch 5k introduced by the host material can 
be compensated through a frequency resonance offset for 
one of the fields. This is in stark contrast with the op- 
timization of OPA for a weak signal in the off-resonant 
case for both ordinary materials and NIMs. Note that 
the requirement of linear phase mismatch is known to 
optimize OPA for strong fields due to the fact that the 
ratio of the intensities of the coupled fields and, conse- 
quently, the coupling phase varies along the medium. A 
strikingly similar behavior occurs for the weak signal and 



homogeneous strong fields in resonant NIMs. To demon- 
strate the outlined extraordinary NLO features and to 
prove the proposed possibilities for producing the tai- 
lored transparency of NIMs and FWM BWMOPO, we 
have adopted the following model for numerical simula- 
tions presented in Fig. [2fc)-(e). The characteristic for 
molecules embedded in a solid host were chosen as: en- 
ergy level relaxation rates r„ = 20, T g = T m = 120; par- 
tial transition probabilities j g i — 7, 7 g „ = 4, j mn = 5, 
jmi = 10 (all in 10 6 s _1 ); homogeneous transition half- 
widths r lg = i, r lm = i.9, r ng = 1.5, r nm = 1.8 (ail 

in 10 12 s- 1 ); T gm = 5, T hl = 1 (all in 10 10 s" 1 ). We 
assumed that A2 = 756 nm and A4 = 480 nm. The 
density-matrix method described in [9J is used for cal- 
culating the intensity-dependent local parameters while 
accounting for the quantum nonlinear interference effects 
(NIE). This allows us to account for changes in absorp- 
tion, amplification, refractive indices and in magnitudes 
and signs of NLO susceptibilities caused by the control 
fields. The changes depend on the population redistri- 
bution over the coupled levels, which strongly depends 
on the ratio of the partial transition probabilities. Fig- 
ure [DJc)-(e) depicts the nonlinear interference structures 
in the local optical parameters entering into Eq. (|3|) at 
the indicated control field strengths and resonance detun- 
ings. Here, Qi = u>\ —u> m l\ other resonance detunings flj 
are defined in a similar way. Coupling Rabi frequencies 
are introduced as G3 = E 3 di g /2h, G 4 = E^d nm j2h. The 
value clnim denotes absorption introduced by the host 
material, and a\o is the fully resonant value of absorption 
introduced by the embedded centers at u>x = uj m i with 
all driving fields turned off. With the given parameters, 
the control fields cause essential changes in the level pop- 
ulations, n « 0.487; r g w 0.484; r n « 0.014; r m « 0.014 
(r; + r„ + r g + r m — 1) which is followed by an appre- 
ciable population inversion at the idle transition. Laser- 
induced NIE in the NLO susceptibilities are followed by 
corresponding changes in absorption and refractive in- 
dices. For the above indicated optical transitions, the 
magnitude G ~ 10 12 s _1 corresponds to control field in- 
tensities of I ~ 10 - 100 kW/(0.1mm) 2 . As discussed 
above, the output signal presents a set of distributed 
feedback-type resonances that provides for BWMOPO. 
Optimization of the output signal at z = is determined 
by the interplay of absorption, idler gain, FWM and by 
the wavevector mismatch. This therefore appears to be 
a multi-parameter problem with sharp resonance depen- 
dencies. The results of numerical analysis of Eq.® for 
one such transmission resonance based on the steady- 
state solutions to the density matrix equations [9j is given 
in Fig.[2^h). The intrinsic absorption of the host slab in 
the NI frequency domain has been assumed equal to 90%. 
In Fig. [2](h), we also introduce the slab length scaled 
to the resonance absorption length, L ra — aj }. The 
value of L/ L ra is proportional to the product of the slab 



length and the number density of the embedded centers. 
Figure [21(f) displays a narrow geometrical resonance at 
L = L r = 32.6427L ra for the optimum frequency offset 
for the signal at Qi/Ti n = y$. This provides the opti- 
mum compensation the given phase mismatch, 8k, intro- 
duced by the host material. Changes in Afe introduced 
by molecules are accounted for within the simulations. 
Fig. EJh) shows that the width of the oscillation reso- 
nance is on the scale of the narrowest (here Raman) tran- 
sition width and the resonant absorption length. Ampli- 
fication in the maxima in Fig. [Hh) reaches many orders 
of magnitude, which indicates the feasibility of cavity- 
less oscillations. Assuming a resonance absorption cross- 
section cio ~ 10 -16 cm 2 , which is typical for dye 
molecules, and a concentration of molecules N ~ 10 19 
cm~ 3 , we estimate aio ~ 10 3 cm -1 and the required slab 
thickness to be in the microscopic scale L ~ (10— 100)/x. 
The contribution to the overall refraction index by the 
impurities is estimated as An < 0.5(A/47r)aio ~ 10~ 3 , 
which essentially does not change the negative refractive 
index. 

In conclusion, we propose mirrorless generation of 
a backward wave and entangled counter-propagating 
left-handed signal and right-handed idler photons in a 
strongly absorbing microscopic samples of negative-index 
metamaterials. The proposal implements the backward- 
ness of electromagnetic waves inherent to NIMs and co- 
herent energy transfer from the control fields to the signal 
through optical parametric coupling. One of the pro- 
posed schemes allows for the independent engineering of 
a negative index and a strong nonlinear-optical response 
in the metamaterial. 
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